Large-scale thermal convection in a horizontal porous layer 
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In a range of physical systems the first instability in Rayleigh-Bernard convection between nearly 
thermally insulating horizontal plates is large-scale. This holds for thermal convection of fluids satu- 
rating porous media. Large-scale thermal convection in a horizontal layer is governed by remarkably 
similar equations both in the presence of a porous matrix and without it, with only one additional 
term for the latter case, which, however, vanishes under certain conditions {e.g., two-dimensional 
flows or infinite Prandtl number). We provide a rigorous derivation of long- wavelength equations 
for a porous layer with inhomogeneous heating and possible pumping. 

PACS numbers: 44.25.-|-f, 44.30.-l-v, 47.54.-r 



Rayleigh-Bernard convection, i.e., a thermal buoyancy 
convection in a layer confined between two rigid horizon- 
tal plates and uniformly heated from below, is one of the 
classical problems of fluid dynamics. When the bounding 
plates are nearly thermally insulating, the first instability 
of the heat conductive state (quiescent state) of the layer 
is long- wavelength (large-scale) [H, 0, H, 01 , «-e., horizontal 
gradients of the fluid velocity field are small compared to 
vertical ones and, thus, the horizontal scale of the flow is 
large compared to the layer height. For large-scale pat- 
terns, temperature perturbations 6 are nearly uniform 
along the vertical coordinate z and in terms of 6{x, y) 
the problem is two-dimensional (2D). Making allowance 
for the heating inhomogeneity, one can find 

dtO+A'^e -V -{Os) +v ■[q{x,y)ve ~ve\ve\'^] =o, (i) 

As = 0, [V xs\, = B=^^[V x{veAe)],. (2) 

Here dt denotes the partial time derivative, q{x, y) is the 
relative departure of the heat flux (imposed by the heat- 
ing) through the boundaries from the critical value for a 
homogeneous heating, Prandtl number Pr = vlx is the 
ratio of the kinematic viscosity and the heat diffusivity, s 
is responsible for advectional heat transfer. Knobloch [Sj 
has discussed pattern selection for Eq. ([T]) with uniform 
q and s = 0, which is relevant, e.g., to Pr ^ oo or 2D 
setups where flows are homogeneous in the y-direction. 
Later on, the term missing in [5|] was accounted for in (|. 
Remarkably, in [7], equations similar to ([1]) have been 
obtained for large-scale turbulence in Rayleigh-Bernard 
convection. 

Various problems related to large-scale thermal con- 
vection in porous media have been repeatedly addressed 
in the literature {e.g., [1, However, researches typ- 
ically either deal with the case of a rectangular cavity 
with a large aspect ratio (0], where no equations similar 
to (|1I2|) have been derived) or discuss a finite departure 
from the stability threshold ([4]; such a departure leads 
to discontinuities of the velocity field, which cannot be 
treated within the framework of the long-wavelength ap- 
proximation). Pattern formation under heating inhomo- 
geneityhas been rather extensively studied for thin films 
{e.g., [8[ and Refs. therein) owing to important technical 



applications (rupture of lubricating films etc.). To 
the authors' knowledge, for porous media the equation 
system similar to (|ll2p has not been presented in the lit- 
erature. 

In this Brief Report, we introduce the specific physical 
system we deal with. Then we derive the equation of 
large-scale convection for this system and find that it is 
similar to Eqs. (|ll2p with B = 0, though the relationships 
between the fiow and the temperature perturbation are 
different for different fluid dynamical systems [1, 0, 0] ■ 

a. Thermal convection in a horizontal porous layer. 
Let us consider convection of a fluid saturating a hori- 
zontal porous layer heated from below. Boundaries are 
impermeable, the heat flux across the layer is fixed (im- 
plying the heat diffusivity of the boundaries is small com- 
pared to the one of the porous matrix saturated with the 
fluid), but inhomogeneous along the layer. Relaxation of 
the local temperature difference between the porous ma- 
trix and the fluid is assumed to be fast, and we do not 
introduce several temperatures for them. For small per- 
turbations of the temperature about Tq, one may guess 
a linear dependence of the fluid density on the temper- 
ature; p{T) ^ po{l - P{T - To)), where po = p{To), 
(3 = {dp/dT)p. The reference frame is such that the 
{x, ?/)-plane is horizontal, z = and z = h are the lower 
and upper layer boundaries, respectively. We adopt the 
conventional Darcy-Boussinesq approximation [2,] , 
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b-^V ■ {vT) 
= 0, d,T = 



+ g(3Te 
-A{1 + 



(3) 
(4) 
(5) 

q{x,y)), (6) 



where H: the macroscopic filtration velocity, K: the per- 
meability, g = —gcz'- the gravity, b: the specific heat 
capacity of the saturated porous medium divided by the 
one of the fluid, x- the heat diffusivity of the saturated 
porous medium, xCpA{l+q{x, yj): the imposed heat flux 
{Cp is the specific heat capacity). 

It is convenient to measure the length by layer height 
h, time by h'^/x^ the velocity by bx/h, the temperature 
by Ah, and the pressure by bp^vxl K . The dimensionless 
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parameter governing the behavior of the system is the 
Rayleigh-Darcy number R — (3Ah? gK /bvx- 
The dimensionless form of system (jSHl]) reads 

-Vp-v + RTe,_^0, (7) 

dtT + V ■ [vT) = AT , (8) 

V • t7 = , (9) 

z-0, 1: v,^0, d,T^^l-q{x,y). (10) 

b. Long-wavelength approximation. For a uniform 
fixed heat flux the first instabihty is known to be long- 
wavelength This holds for q{x,y) slowly varying in 
the boundary plane: |Vg|/|g| ~ e <C 1. In order to avoid 
large temperature gradients, which correspond to jumps 
of derivatives of fields in the long- wavelength limit 4], 
we restrict ourselves to a small supercriticality. Below we 
will find the critical value Rq for a homogeneous heating 
and set R = Rc', hence, nearcritical regimes occur for 
small q{x,y) which may be thought to measure the rel- 
ative departure of the Rayleigh-Darcy number from the 
critical value. 

Owing to Eq. ([9]) the horizontal component of velocity 
V is large compared to the vertical one because variations 
of velocity H by shifts along the horizontal directions are 
small compared to ones by shifts transversal to the layer. 
Let us explicitly account for this fact, v = wez + s~^u, 
where m is the horizontal component of the filtration ve- 
locity field, w is the vertical one. Rescaling horizon- 
tal coordinates, (x^y) {£~^x, £~^y), explicitly writing 
g(x, y) — £^q2{x, y), and projecting momentum conserva- 
tion law ([7]) onto the vertical and horizontal directions, 
one may rewrite equation system ((Tl llOp in a form con- 
venient for the further treatment; 

~dzp-w + RT = 0, (11) 

u^-e^^2P, (12) 

dtT + d,{wT) + V2 • [uT) = dlT H- £2A2T , (13) 

d,w^V2-u = 0, (14) 

2 = 0,1: w^O, d,T^^l~£^q2{x,y). [lb] 

Here the subscript "2" for spatial derivatives means the 
differentiation with respect to two horizontal coordinates. 

Since e appears squared in piT[T5|) . only even powers 
of e are present in the expansion: w = 'Wo + e'^'W2+£'^W4^ + 
. . . , r = To + e'^T2 + e^T4 + . . . , etc. The long- wavelength 
approximation assumes a weak spatial inhomogeneity of 
temperature perturbations which results in a slow tem- 
poral evolution. With only even powers of e in the ex- 
pansion, one expects characteristic times of the evolution 
of large-scale patterns to be not lesser than oc there- 
fore, in terms of "slow" times, dt = e^^ta + £^f^t4 + • • ■ ■ 
£° : In the leading order, problem (fTTHTSl) yields 



-d,po - Wo + RTo ^ , (16) 

uo^O, (17) 

a. (woTo) + V2 • {uoTo) = din , (18) 

d,wo -h V2 • wo - , (19) 

z = 0,l : wo = 0, d,To = -1 . (20) 



From Eqs. (|19ll7p . d^WQ = 0, = Ci{x, y) = 
[Ci — due to the boundary conditions (BCs)]. 
From Eq.lUHl), d^To = 0, i.e., Tq = C2{x,y)z + e{x,y); 
accounting for BCs (|20p . one obtains 

To = -z + 0ix,y), (21) 

where 9{x,y) should be determined from higher orders 
of the expansion while here it appears as an unknown 
fimction of the horizontal coordinates. 
From Eq. mD, d^po ^ RTq ^ -Rz + Re{x,y), 

Po^-^Rz'' + R0ix,y)z + noix,y), (22) 

where IIq^x, y) is unknown in this order of the expansion, 
e!: -d,p2-W2+RT2 = 0, (23) 

U2 - , (24) 

dt.To + d,{w2To) + V2-(u2To) = dlT2 + A2T0 , (25) 
d,W2 +\72-u2=0, (26) 
z = 0, 1: W2=0, d,T2 = -q2{x,y). (27) 

From Eq. (IM)) . 

U2 = -V2P0 = -RzV29{x,y) - V2no(x,y). 

From Eq. 

dzW2 = -V2 • U2 ^ RzA29{x,y) + A2i7o(a;,y), 
W2 = ^Rz^A29{x,y) + A2na{x,y)z + C3{x,y). 

BCs (HZl) for the velocity field yield 

z = : C3 = , 

z = l: A2no{x,y) = ^^RA29{x,y). 
From the latter BC, 
no{x,y) = ~\R9{x,y) + Tro{x,y) , A27ro(a::,y) = . 

Note, {U2) = — V27ro (henceforth, (/) = fdz). Let 
us consider a layer domain limited in the horizontal di- 
rections by boundary F (to keep it simple, we assume 
r to be vertical), where a fixed fluid gross flux (or its 
absence; "gross" means averaged over z, and, in partic- 
ular, the absence of the gross flux does not necessarily 
claim the absence of the flux) is imposed. Mathemat- 
ically, this means that the orthogonal to F component 
of (-1*2) is flxed: (w2)„ = Q (Q is not uniform along F; 
owing to mass conservation, Jj^QdF = 0). Hence, one 
obtains the boundary problem for no(x,y), 

A2^o = 0, dnno\r^-Q (/^Qdr = 0), (28) 

where (?„ is the orthogonal to F component of the gra- 
dient. This problem has a unique solution (up to an in- 
significant constant) unambiguously defined by Q. Thus, 
ttq{x, y) describes an imposed advection (pumping) in the 
layer, which is caused and unambiguously controlled by 
BCs on F (the pressure or the fixed gross flux). As soon 
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as we consider a nearcritical behavior, it makes sense to 
not allow for an imposed advection in the leading order, 
otherwise a moderately strong advection would overpress 
the effect of a weakly inhomogeneous heating. We set 

TTo{x,y)^0 (29) 

and will take the imposed advection into account in 
higher orders. 

With (|29p . one can write down 

U2 ^-R{z- i) V20(x, y) , (30) 
W2^^Riz^-z)A29ix,y). (31) 

Let us now integrate (pS)) over z E [0, 1]; 

dtATo) + {d,{w2To)) + V2-{u2To) = (a^Ta) + A2 (To) . 

(32) 

Due to BCs (HZl), J^d,{w2To)dz = W2To\l^o = 0- 
a^Tadz = d,T2\l^o - -q2{x,y)\UQ = 0. Whereas, 

V2 • {U2T0) - -V2 ■ {U2Z) + V2 • {{U2) e) . 

Substituting ([30]) into the first term of the rhs of the last 
equation and accounting for {U2) = — V27ro = 0, one finds 

V2 • (W2T0) = RA20{z^ - z/2) = (i?/12)A20. 

For the rest of the terms in Eg. (|32|) . dt^To = dt^O, 
A2T0 = A26', and Eq. ^ finally reads 

dt^e = {I -R/ 12)^29. (33) 

For R < 12, Eq. (|33p is a conventional diffusion equa- 
tion, and for trivial BCs (on F) or an infinite layer, all 
inhomogeneities of 6 decay. For i? > 12, it is a diffu- 
sion equation with a negative diffusivity, where all the 
inhomogeneous perturbations grow. Thus, i? = 12 is the 
linear stability threshold of the system. Note, nonlinear- 
ity does not play a role in this order of the expansion. 
In order to account for nonlinear effects and the depen- 
dence of the linear stability on the wavelength (now all 
the perturbations either grow or decay regardless to their 
wavelength), we should restrict ourselves to the vicinity 
of the stability threshold. For this purpose, we set 

R = Rc = 12 

and introduce departure from the threshold via q. Since 
the local Rayleigh-Darcy number ( "local" means defined 
for a small domain of the layer) i?iocai = R{^+q), positive 
q corresponds to a supercritical regime, negative q does to 
a subcritical one. As soon as R= 12, Stjfl = 0; therefore, 
we should consider a slower evolution, dtO = e'^dt^O + . . . . 
Let us now derive T2 from Eq. ([^5]) . 

dlT2 = d,{w2To) + V2-(m2To) - A2T0 ; 
T2 = - (|z4 - 2z3) A26I + (2z3 _ 3z2)6'A26' 

+ {z^ ~ z^)A29 - (2z3 - 3z^)\/2-{0V2d) - \z'^A2B 
+C4{x,y) z + e2{x,y) . 



Due to the relation V2 • {0V29) = 9 A29 + {^20)'^ , one can 
eliminate the term 9A29 from the expression for T2; 

T2= (-iz4 + ^3_ 1^2) 

+ (-2z3 + 3z2)(V20)2 + C4{x, y)z + 92{x,y). (34) 

Owing to BCs (gT]), Ci{x,y) -172(3;, y); 6*2(2;, y) is still 
undefined. 

Remarkably, 9 and 02 depend on z in one and the same 
fashion (i.e., are uniform along z) and are for the mo- 
ment undetermined functions of x and y. Hence, 02 may 
be chosen as one needs, and this will be automatically 
balanced by 0. Let us use this fact. Notice, 

l^T) = -l + 0{x,y) + e^ {T2) + 0{e^)- 

therefore, if one defines 02 in such a way as to make 
(T2) = 0, then will be a z-mean temperature up to the 
truncation accuracy of our expansion, i.e., e**. Thus, 

{T2) = -^A20 + i(V20)2 - ig2(a;, y) + 92{x, y) = , 

62 ^ ^A29 ~ \{V20f + \q2 ■ (35) 

Finally, 

T2^{-\z^ + z^-\z^ + h)^^0 

+ (-2z3 + 3z2 - i) {V20f + {-z + \) q2. (36) 

From Eq. follows dzP2 = -W2 + RT2, and, inte- 
grating it with respect to z, one finds the pressure 

P2 = (-§2^ + 3z4 - 4z3 + 3z2 + iz) A29 
+ (-6z4 + 12z3 - 6z) {V20Y 
+ (_6z2 + 6z)g2 + il2(x,y). (37) 

ef : As we will not construct the expansion beyond this 
order, we do not have to consider all the equations. The 
following is sufficient; 

M4 -V2J32 , (38) 
duTa+d,{w4To) + V2-{uiTo) 

+dAw2T2) + V2-(u2T2) = d^.Ti + A2T2 , (39) 
d,wi + V2 ■ U4 = , (40) 
z = 0, 1 : W4 = 0, d^Ti = . (41) 

Without calculating M4, we may substitute (|38p di- 
rectly into (gni); dzW4 — A2P2- Then 

y04 = (-iz6 -I- |z5 -z^ + z^ + ±z^) Ap 

+ (-fz5 + 3z4-3z2) A2(V2e)2 

+ (-2z3 + 3z2) A2q2 + ZA27T2 . (42) 

From BCsgl]), 

wil^i = - \A10 - lA2{V20f + A2g2 + A27T2 ; 

therefore, 

n2^^\A20+l{V20? -q2+TT2, 
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where A27r2 = 0. 

Notice, making use of Eq. ([57]) , one can calculate 

("4) = -V2 (P2) = -V27r2(a;,?/) . 

We have already established the relationship between 
{U2) and ttq and found ttq to describe an imposed ad- 
vection. We have set tto = 0, but now it makes sense 
to keep 7r2, since it describes an advection imposed by 
lateral boundary conditions. 

Substituting 772 into Eq. (|42p . one obtains the final ex- 
pression 

= (-1^6 + 3^5 _ ^4 ^ ^3 ^ ^^2 _ 1^) 
+ (-§ z5 + _ 3^2 ^ A2(V2(?)2 

+ (-2z3 + 3z2 - z) A2g2 . (43) 

Now the integration of Eq. p9|) over z £ [0, 1] yields 
the evolution equation for 9 

+ {d.{w2T2)) + V2-(W2T2) = {dlTi) + A2 (T2) . 

The mean values of all z-derivatives are zero due to BCs; 
additionally, (T2) = 0. The remainder is 

due + V2- {UATQ + U2T2) ^0. (44) 

With To, T2, P2, and U2 known, one can find 

V2 • (uiTo) = -V2 • (To V2P2) = ^ A20-V2 • (0V27r2) ; (45) 

V2-{U2T2) - -V2-(|V20(V20)2 - g2V20) . (46) 

Substituting (|45l) and ([461) into Eq. (gH) , one obtains 
the slow evolution equation for in the final form which 
is similar to Eq. ([1]); 

dtO + U-W9+ ^A^e - V- (I V6i(V6i)2 - qV0) = 0. (47) 

Here the imposed advection U = — Vvr — (u), \7 -U — 0, 
subscripts "2" for differential operators are omitted as 



all the fields depend on x and y only, indexes indicating 
the smallness order are also omitted as this equation re- 
mains valid in original terms (without formal smallness 
parameter e). From Eqs. (|30l3f p the fluid flow up to the 
leading order of accuracy is 

w = 6(1 - 2z)V6'(a;, y) + 6{z^ - z) Ae{x, y)e^ . 

Advection speed ?7 cx e'^ is small against u oc e when 
flow is excited, but is important due to its properties: 
in contrast to u, U provides a nonzero gross fluid flux 
through a vertical cross-section of the layer. 

Note, the heating inhomogeneity makes the quiescent 
state impossible, and from Eqs. (l36|) and (|43|) it follows 
that below the convective instability threshold, when 9 
decays to zero, the establishing state has nontrivial 



«;bg = (-2z3+3z2-z)Ag, Tbg = {^-z) g - z. (48) 
However, these fields are small against the fields excited 
above the threshold. 



Remarkably, Eq. (|47p holds for the case of a uniform 
heating and a weakly inhomogcneous porous matrix. In 
the most general case, q{x, y) should be replaced with 



qgen{x,y) = q{x,y) 



K{x,y) boxa 
b{x,y)xix,y) Kq 
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[meanwhile, the background state (|48| for this general 
case cannot be obtained via the plain replacement of q 
by qgcn]- This remark is important for experiments as it 
is more convenient to maintain/control a uniform heating 
of the layer with a weak inhomogeneity of the structure 
of the porous matrix, which is inevitable in real systems. 
Owing to the same reason, the physical system we have 
discussed is the most relevant one for works address- 
ing phenomena related to spatially localized convective 
currents excited under parametric disorder [frozen ran- 
domly inhomogcneous q — q{x)]. 
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